function f=find_math_bi(h,kk1,Y_i)

if kk1==1
	f=2.^(-1/2).*h;
else
	f=(-15).*h.^(-4).*(kk1-1).^(-5).*pi.^(-5).*cos((kk1-1).*pi.*Y_i).*(3.*h.*(kk1-1).*pi.* ...
		cos(h.*(kk1-1).*pi)+((-3)+h.^2.*(kk1-1).^2.*pi.^2).*sin(h.*(kk1-1).*pi));	
end

	%
